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Limits on the size of quenched chiral logs in the pion mass for Wilson fermions are investigated. The smallness 
of chiral logs is shown to be a result of the suppression of the hairpin diagram for small p^, such that the value of 
the hairpin on the pion mass shell is much smaller than the physical m^, . A direct calculation of the topological 
susceptibility from the same data gives m^/ « 1 GeV. 



1. INTRODUCTION 

One of the most direct manifestations of closed 
quark loop effects is in the mass of the flavor sin- 
glet pseudoscalar rj' meson. The rj' propagator in- 
cludes not only the valence quark-antiquark term 
that appears in the nonsinglet propagator, but 
also the contribution from diagrams in which the 
valence quark and antiquark annihilate. These 
"hairpin" diagrams probe the topological struc- 
ture of the gauge field via the axial U{1) anomaly. 
In a large Nc approximation, the axial anomaly 
may be introduced perturbatively as a term which 
breaks the U{1)a symmetry of the massless quark 
Lagrangian. From the chiral Lagrangian point 
of view, the anomaly adds a term Ci to the 
U{Nf) X U{Nf) Lagrangian Co which gives the rj' 
its mass. For low momentum, and keeping only 
terms quadratic in the rj' field, the most general 
form of such a term is (c.L Ref. 0) 

^i = \ {{A 1) [{d.rjT - mlv")] - Aml^'')il) 



This combines with the term in Co to give an rj' 
mass of rrir]' — \/ ttiq + mj. The constant A is 
a renormalization of the rj' field induced by the 
inclusion of the anomaly. In the large Nc frame- 
work, the term (|l|) may be identified with the 
quenched hairpin diagram. It corresponds to an 
amputated hairpin vertex of the form 



The main difference between quenched and full 
QCD is that in quenched QCD, the hairpin ver- 
tex appears only once in the t]' propagator, while 
in full QCD it appears an arbitrary number of 
times. In the latter case, the hairpin insertions 
sum up geometrically and shift the = Gold- 



stone pole to a pole at p 
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By contrast. 



the quenched rj' propagator includes only a sin- 
gle hairpin insertion. Not only is the Goldstone 
pole not cancelled, but the hairpin graph adds a 
double pole l/(p^ — wj)^ term to the propagator. 
Note that, by (|^), the quenched diagram with a 
single hairpin insertion includes both a single pole 
and a double pole term, with coefficients 1—A and 
ArriQ respectively. 

The appearance of a double pole in the 
quenched rj' propagator gives rise to anomalous 
chiral behavior, e.g. in the relation between the 
pion mass and the quark mass[0j|]. The chiral 
symmetry result that is linear in the quark 
mass is replaced in the quenched approximation 
by 
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where the parameter 6 which determines the 
anomalous power behavior is the coefficient of the 
quenched chiral log in the one-loop graph, and is 
proportional to Am^, the value of the hairpin in- 



sertion at p^ w mj. This gives 



Ii{p'') = -{A-l){p'-ml) + Aml 



(2) 



At 



24^2/2 



(4) 



2 



If we assume that ^4 w 1, a rough estimate using 
mo ~ 0.9 GeV gives 

6^0.2 (5) 

2. QUENCHED LOGS IN THE PION 
MASS 

It has recently been argued that the behav- 
ior of the pion mass as a function of the bare 
quark mass calculated in quenched lattice QCD 
shows little or no evidence for the presence of chi- 
ral logs at the level suggested by the estimate (|^) . 
We have analyzed the quenched pion mass from 
ACPMAPS data at four different /3 values and a 
variety of hopping parameters, as shown in Ta- 
ble 1. The pion masses were extracted from the 
pseudoscalar propagator with pointlike sources. 
To determine and remove the effect of excited 
states, both one- and two-exponential fits for a 
variety of time windows were carried out. For 
some of the k values at (3 = 5.7, the pion mass 
obtained in this way was compared with that ob- 
tained using smeared-source quark propagators, 
and the results were found to agree within sta- 
tistical errors. For each /3 value, the pion masses 
for N values of k were calculated (here iV = 3 or 
4), with the full N x N error matrix being com- 
puted by a jackknife elimination. By minimizing 
the covariant x^, a 3-parameter fit was obtained 
to the fitting function 

ml ^ C {k-' - ^-')^ (6) 

with fit parameters C, Kc, and 5. The results for 
the parameter S are given in Table 1. For all four 
values of f3, the value obtained for S is consistent 
with zero. Combining the statistics of the four 
/3's, we get 

(5 = 0.00 ±.03 (7) 

In the analysis leading to (0), only values of k cor- 
responding to pion masses of less than 750 MeV 
were included in the fits in an effort to minimize 
the effect of higher order chiral perturbation the- 
ory terms. In the f3 — 5.7 result (where there 
were four mass points and hence one degree of 
freedom in the fit), the covariant was 0.1, in- 
dicating that a good fit is obtained without the 
need for higher order terms in (ph. 



Table 1 



Exponent S from vs. niq. 
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k's 


6 


5.7 


.161, .165, .1667, .168 


.015(47) 


5.9 


.157, .158, .159 


-.004(62) 


6.1 


.153, .154, .1545 


.009(75) 


6.3 


.1510,. 1513,. 1515 


-.033(114) 



3. CALCULATION OF THE HAIRPIN 
DIAGRAM 

To investigate the apparent suppression of chi- 
ral logs further, we calculated the hairpin dia- 
gram directly at /3 = 5.7, using the technique of 
Kuramashi, et al Q. The calculations were car- 
ried out on both 12^ x 24 and 16^ x 32 lattices. 
Our main conclusion is that the value of Ara^ 
extracted from a direct calculation of the hairpin 
graph is much smaller than the physical rj' mass- 
squared, and consistent with the small value of 5 
inferred from the limits on quenched chiral logs in 
the pion mass. At the n values and lattice sizes 
for which a direct comparison could be made, our 
raw data was in good agreement with that of Ref. 
Q . Our somewhat different conclusion regarding 
the suppression of the hairpin vertex at small p^ 
follows from several factors which we briefly men- 
tion here. First, in order to extract the coefficient 
of the double pole term in the hairpin (i.e. Atoq), 
the pion propagators on either side of the vertex 
were assumed to contain excited-state as well as 
ground-state contributions, with the excited state 
mass determined from the previously described 
pion propagator analysis. (The coefficient of the 
excited state term was a free parameter in the 
fit.) It is important to note that, in the hairpin 
diagram, terms involving an excited state on one 
but not both sides of the hairpin still contain a 
single Goldstone pole, and thus fall off with the 
same exponential factor as the ground state term. 
They are only suppressed by a power of t, i.e. 
they fall off like e"™"* instead of ie~™"*. For our 
fits, the inclusion of excited state contributions 
to the hairpin propagator reduced the extracted 
value of Atuq by about 30 to 50% compared to 
a pure ground state fit. Secondly, a comparison 
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♦ p = 5.7 12^x24 
o 16^x32 



inula, which relates uiq to Xt , the topological sus- 
ceptibility of pure glue, 
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Figure 1. The value of the chiral log parameter S 
extracted from the hairpin calculation (using S = 
Auiq / 247r'^ f^) compared with the one standard 
deviation upper bound from vs. niq. 



of our results for 12'^ x 24 and 16"^ x 32 lattices 
with the same /3 and k values exhibits a large 
and highly mass dependent finite volume effect 
on the measured value of ArriQ. This raises doubt 
about the results of performing a chiral extrap- 
olation on a fixed size box. (The chiral extrap- 
olation contributed substantially to the quoted 
value of mo « 750 MeV in Ref. g.) The re- 
sults of the hairpin calculation for several k's and 
two box sizes are converted to an effective value 
of S and plotted in Fig. 1. They are seen to be 
consistent with the one-standard-deviation upper 
bound from the pion mass analysis. (The dotted 
line represents the upper bound from the /? = 5.7 
pion mass data in Table I.) 

For the range of quark masses considered, the 
results presented here provide strong evidence 
that the size of quenched chiral logs is suppressed 
by the fact that the hairpin vertex evaluated on 
the pion mass shell is much smaller than expected 
from the assumptions that uiq « rn^, and A ^ 1. 
If we discard the possibility that niQ << m"^, 
(which would be a disturbing failure of QCD to 
reproduce the real world), it may be concluded 
that A « 1 and that the hairpin is highly mo- 
mentum dependent. In addition to analyzing the 
time-dependence of the hairpin propagator, there 
is an indirect way to determine mg independently 
of A by appealing to the Witten-Veneziano for- 
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It turns out that the same data generated in the 
hairpin calculation can also be used to obtain an 
approximate measurement of the winding number 
of each gauge configuration in an ensemble, using 
the anomalous chiral Ward identity. 



d x{^Pj^iP)g — — 
m„ 



(9) 



By studying the behavior of a single 7^ loop in- 
tegrated over the entire lattice as a function of 
quark mass, (^ may be used to obtain an ap- 
proximate determination of the winding numbers 
of the configurations. (A more detailed discus- 
sion of this method and its applications will be 
presented elsewhere.) From the same 7^ loops 
used to compute the hairpin at f3 — 5.7 and 
V = 12^ X 2Aa^, we obtain a mean squared wind- 
ing number of {v^) — 23 ± 3. The error here is 
a very rough estimate based on varying the cri- 
teria for observing a l/ruq pole. From this, and 
using a^^ = 1.15 GeV, we obtain the topological 
susceptibility 

Xt = {i^^)/V « (ISOMeV)^ (10) 

and from the formula (||), toq ~ 1.1 ± 0.2 GeV. 
This is consistent with results obtained by the 
cooling method 

We are grateful to George Hockney for many 
contributions to this research. 
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